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We deduce the particle distributions in particle collisions with 
multihadron-production in the framework of mechanical statistics. They 
are derived as functions of x, P^ and the rest mass of different species for a 



a.: 

^ (— i fixed total number of all produced particles, inelasticity and total transverse 



energy. 



H ' For Pt larger than the mass of each particle we get the behaviour 

1 ~ yJPre t h (1) 



dP T 

Values of < Pt >tt, < Pt >k, and < Pt >p in agreement with experiment 
are found by taking Tjj = 180MeV (the Hagedorn temperature). 

PACS: 13.85.Hd - Inelastic scattering: many-particle final states 
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Predictions of multihadron-production beyond accelerator measuraments is an im- 
portant tool for study of cosmic rays with high energies. Thermodynamical concepts 
have been advocated for the study of these processes up to Tevatron energies [1]. A non- 
isotropic phase space should be considered, since most of the energy of the final state 
comes from the longitudinal momentum of the final particles. This has motivated us in 
a previous work [2] to look for a statistical interpretation of the semi-inclusive statistical 
model developed by one of us [3], by approximating the energy of the final particles 



Eq.(2) follows from the assumption for a dominance of the longitudinal contribution to 
the total energy and is well verified for ultra-relativistic particles at small angles with 
respect to the direction of the initial particles. 

We assume that the missing part of the total longitudinal momentum is responsible 
for creation of particles and giving them a transverse momentum. 

Here we perform a mechanical statistical description of this picture by solving the 
problem of partition of energy in N produced hadrons, sharing it in the final longitudinal 
momentum and the total transverse energy. 

By the classical method of the statistical mechanics one finds for the distribution of 
the various particles (in each hemi-sphere): 



where k and the dimensional constants M 2 and /Xq are determinated by the boundary 
conditions: 






(3) 




(4) 




(5) 




(6) 
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Here the sums exclude the leading particle and r] ine i and E T are the inelasticity and the 
total transverse energy, respectively (the latter has a negligible contribution with respect 
to the longitudinal one in Eq.(5)). 

From Eq.(3) it is easy to derive the total number of particles of the species i: 
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rii = J dxe kx e (7) 

by keeping into account that the function to be integrated has a sharp maximum around 
the minimum, x™ m = m*/ )j,oVk, of 

2 

f l (x)=kx + ^ (8) 

X/Iq 

and fi(x) may be well approximated by 

u{xT m ) + \r:{xT m ){x-xT m ) 2 (9) 

which gives rise to a gaussian integral. 

In Table I we report fi(x™ m ) and f"{x™ %n ) and the corresponding values for 

the functions 

9 f(x) = f t (x)-p\n(x) (10) 

which appear in the expression of < >, < P^ >= A*o < x i > an( ^ °^ ^ ne contribution of 
the term proportional to mf to the total transverse energy in correspondence of p — 1/2, 
1 and — 1, respectively. 

From Table I it is easy to derive within the gaussian approximation for the integrals 
in the variable x: 



H 



T 4 / / m 2 \ e > 

n t < Pi >= mul < x t >= 4v^^ I 1 + V 1 + # ) 7 ^TVi (13) 



n t < >= (14) 
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ni < _^ >= ^2^-^=- (15) 



where 



The expression in eq.(3) gives the correlation known gul" effect (with the square 

of the transverse momentum increasing with x). It also implies a larger transverse mo- 
mentum for the heavier particles, which are produced at a larger x due to the presence of 
the exponential factor exp^—mf/xn^). 
The ratio 

T7sr"» < 17 > 

does not depend on the species and this property is an important check of this statistical 
approach. By integrating eq.(3) in x, again with the gaussian approximation for the 
function, one obtains from g _1 (:r) by the substitution m 2 — > m 2 + Pj, , one gets the Pt 
distributions: 



drii 2\/2txPt 1 e 



1 + 



rp'Z 

1 H 



(18) 



which, for Pt { ^> mj, T# approaches: 



dn t 2y/2nP Ti T H _ P J± 

= — - e H 

dP T . M 2 
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(19) 



So we recover the exponential behaviour found experimentally. By fixing T H = 180MeV 
(the Hagedorn temperature [4]), we predict < Pt z > as a function of rrii and T H : 



3/2 A 



K ( I 4- m A 



which gives: 



< P T > 7r = 390MeV (21) 

< p T >K = AQ7MeV (22) 

< p T >_= 580MeV (23) 

to be compared with the experimental results at s = (5A0GeV) 2 , 372MeV, A82MeV and 
QOQMeV respectively [5]. Also we find: 



-^ = J— e t h = 0.256 24 



rip rriN 



^ = 2j—e t h = 0.061 (25) 

while the experimental values are 0.112 ± 0.010 and 0.078 ± 0.011, respectively [5]. The 
approximation described in eq.(2) is valid as long the P T distribution is considered for: 

(/ m 2 + P 2 \ 
yi+ l T 2 R Ti - 1 (26) 

which is very well satisfied for any value of rrii and P T ., if one takes k = 16.5, as it follows 
from the value 7 found for k by [3] at s = (20GeV) 2 by assuming: 

k(s) = k(s ) (£f 3 (27) 

The high value of the second derivatives in Table 1 support the validity of the gaussian 
approximation. Also the functions fi(x) and g ip (x) are not so different, since they differ 



by a term proporional to ln(rr), and so we expect a better precision for the ratios of the 
quantities described in eqs. (11-13), which give < Pr t >, < > and < Xi >. From 
eq.(16) with T H = 180GeV and k = 16.5, one finds at ^ = 540Ge1/, /i ~ lA6GeV, 
which implies for the average contribution to the transverse energy: 

< >= ^ ~ AMeV (28) 

Xy/S V s 

for the transverse momentum (independent of the species) and 



2 2 m,- + -»/T^+m? 

which gives 0.56MCV, 4AMeV and 18AMeV for 7r, X and AT, respectively. 

The good results found for Pt distributions support the statistical approach de- 
scribed here. The large value found for nx/n^ and slightly small for np/n w may be 
related to have not considered the fluctuations in the r.h.s.'s of eqs. (4-6). In particular 
we know that the number of particles produced at a given energy is not even a gaussian 
(or a Poissonian), as realized from the distribution of the relative multiplicity n/ < n > 
found many years ago [3] and being also analysed in [6]. In fact we may expect that the 
tail at large n will favour a copious pion production, while for small n antiprotons will 
be less disfavoured by the factor e -m i/ x/1 ° by the larger x involved. For the purpose of 
cosmic ray study we have to find the energy dependence of the quantities we have fixed 
(n, Vinei and Et) or of the associated Lagrangian multipliers (k, //) and M 2 , the dimen- 
sional constant related to n. For n we may assume the power behaviour (s/so) a with 
a = 0.13 succesfully assumed in the semi-inclusive statistical model [3] and confirmed by 
the empirical expression found at Tevatron experiments (= —7 + 7.2(s/s )°' 127±0 ' 005 ) [7]. 
For k we could assume, as we have already done in checking the validity of eq.(2), the 
same behaviour implied by the use of the variable x s = x(s/so) a [8]. Finally the energy 
dependence of T# = may be found by studying the high Pt behaviour of dfii/dPr 
as given by eq.(20), as a function of energy: the increase with s of the mean transverse 
momenta implies a corresponding behaviour [5] of T H according to eq.(20). A realistic 
description will require the consideration of the fluctuations. 
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The success of the statistical approach developed for the production of strange par- 
ticles in e + e~ reactions [9] makes us confident in applying the approach described here to 
the production of the particles composed of the light SU(3) quarks (u, d, s). However 
the event found in the associated charm production in deep-inelastic neutrino induced 
reactions at CHORUS [10] may be a signal for a larger production than the evaluations 
based on elementary QCD processes with amplitudes proportional to 0:5 and motivate 
the extension to the charmed particles of the statistical approach described here. 
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Table I 
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